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Characterization of phase transition in Heisenberg mixtures from density functional 
theory 


L. S. Li * and X. S. Chen İ 
Institute of Theoretical Physics, Chinese Academy of Sciences, P.O. Box 2785, Beijing 100080, China 


The phase transition of hard-sphere Heisenberg and Neutral Hard spheres mixture fluids has been 
investigated with the density functional theory in mean-field approximation (MF). The matrix of 
second derivatives of the grand canonical potential Q with respect to the total density, concentration, 
and the magnetization fluctuations has been investigated and diagonalized. The zero of the smallest 
eigenvalue As signalizes the phase instability and the related eigenvector x; characterizes this phase 
transition. We find a Curie line where the order parameter is pure magnetization and a mixed 
spinodal where the order parameter is a mixture of total density, concentration, and magnetization. 
Although in the fixed total number density or temperature sections the obtained spinodal diagrams 
are quite similar topology, the predominant phase instabilities are considerable different by analyzing 
x, in density-concentration-magnetization fluctuations space. Furthermore the spinodal diagrams 
in the different fixed concentration are topologically different. 


PACS numbers: 05.70.Jk, 61.20.Gy, 64.60.Cn, 61.25.Em 


I. INTRODUCTION 


The study of the properties of the Heisenberg liquid 
emerged as one of the fascinating theoretical problems 
since it was stimulated by MC simulations of Heisenberg 
fluid E, which strongly hint that the ferromagnetic tran- 
sition of the Heisenberg fluid challenge the traditional 
viewpoint to ferromagnetic transition theory. Nijmei- 
jer et al. found the value of critical exponents of the 
Heisenberg and Ising fluid |] differ from the expected re- 
sults for the lattice models. They suggested that the the 
effective exponents is related to Fisher renormalization 
with the fixed density. Recently Mryglod et al. sim- 
ulated the larger Heisenberg fluid system and obtained 
the critical exponents from standard finite size scaling 
theory [3] They argued the effective exponents is cor- 
rected by Fisher renormalization when a thermdynam- 
ics system under a constraint [4]. Although the spin liq- 
uids exhibit interesting phase transition by simulations, 
the formal theory of critical phenomena has not solved 
this problem completely. In fact the spin fluids have a 
complex phase behavior as the coupling between the ad- 
ditional spin degrees and spatial coordinates. Besides 
ordinary gas-liquid phase, such a spin fluid model dis- 
plays paramagnetic gas-ferromagnetic Liquid phase tran- 
sitions, critical end point, and tri-critical point. Many 
models have been proposed such as the discrete Ising 
A IG continuous xep 10] and Heisenberg fluid 
, L LA (13) (14) EEI 16, 17]. Tavares et al. us- 
ing both a mean field (MF) and a more refined modified 
mean field (MMF) density functional theory have found 
in some regime mixed first-order transition, namely a 
condensation-ordering transition[12 in the pure Heisen- 
berg fluid. Both theoretical works and MC results showed 
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that the first-order transition (i.e., an isotropic vapor 
phase and a ferromagnetic liquid phase) are the mixed 
transition of the ordering fluctuations and density fluc- 
tuations and Li et al. investigated it using the method 
proposed by Chen et al. is, to characterize the mixed 
phase transition[17]. 

The phase diagrams of binary spin mixture calculated 
by the mean field theory and Monte Carlo simulation 
show many fascinating phase behaviors which come from 
competition among magnetic, condensation and concen- 
tration fluctuations|20}. So far the microscopically moti- 
vated studies try to comprehend the picture of phase be- 
havior of these fluids. However, as we shall demonstrate 
here, knowledge of the complete phase diagram is essen- 
tial but not enough. In order to understand the existence 
of the mixed transition, we have to classify the phase 
transition characters. So our purpose in the present work 
is to characterize the phase transition along ordering- 
condensation-demixing transition line. Therefore, we in- 
vestigate a mixture of Heisenberg hard spheres and Neu- 
tral hard spheres (HHNH) using the method|[18, to 
characterize phase transition. We address these ques- 
tions using density functional theory in the so-called sim- 
ple mean-field theory. As a result, the phase transitions 
of different character take place in this system, as stud- 
ies of the phase diagram and of the associated critical 
phenomena have been shown. 


II. MODEL 


In this paper we consider a binary mixture fluid of 
Heisenberg Hard-Sphere (HHS) and Neutral Hard-Sphere 
(NHS) with equal diameters ø. The pair potential for two 
of such particles at positions rı and rz is given by 


Uablr1, W1, r2, w2) = A(ab)u**(r, w1, w2) + u(r), (1) 


where r = |r] — rə| is the distance between particles at 
positions rı and rə and w; = (6;, Qi) are the orientations 
of particles with 0 < 0; < m and 0 < ¢; < 27. The sub- 
scripts a and b denote the components in fluid considered 
{a(b)= HHS, NHS}, and we define 


(as) ={ 9 a =b = HHS, (2) 
0, a=b=NHS, afb, 

which means all the particles have a diameter yielding 
a repulsive hard sphere potential u” and only the HHS 
particles attract each other via a Heisenberg potential 
uS. The repulsive hard sphere interaction and the spin 
part are given by 


+oo, r< 
hs , a0, 
= 3 
w { 0, r>o. (3) 
and 
0, r<o, 
ut (r,wi, w2) = 4 —J(r)si-s2, o<r<re (4) 
0, Tr>Te 
where 
—z(r/o—1) 


here s; is a unit vector (|s;| = 1) in the direction of the 
spin moment with 


S1'S2 = Cos(w1, w2) = cos 6; cos O2+sin 9; sin 62 cos(ġ1ı—¢2). 


(6) 

In this paper we set the cutoff distance of the Yukawa 

potential re for re = co and the dimensionless parameter 

z for z = 1. To have ferromagnetic phase favoring parallel 

orientations the coupling constant e here is taken to be 
positive. 


Ill. DENSITY FUNCTIONAL AND 
MEAN-FIELD APPROXIMATION 


The grand potential free energy 2 of a nonuniform spin 
liquid mixtures is the minimum of the functional, 


(r), V, T] = Flou(r,w), ou) (7) 
é J dusdep (0, w) (VE, (£,w) — 


4 J dron (r)(VA; (x) 


Olen (r, w), PN 


HH) 


— un), 


where index H and N means HHS and NHS, respectively. 
H is the chemical potential and the intrinsic Helmholtz 
free energy of the inhomogeneous fluid F [py (r,w), pn (r)] 


is a unique functional of the densities py(r,w) and py (r). 
The Helmholtz free energy can be written as 


Flpn(t,~), pu(r)] = F" lon, w) pu()] (8) 
1 


1 
+ >| a f drıdroədwıdw2g(£1, w1, r2, w2; À) 
0 


pH(¥1,W1)u*(r, w1, w2)PH (r2, w2) 

where F*s[pop(r,w),pn(r)] is the Helmholtz free en- 
ergy of Hard-Sphere (HS) mixtures system and 
g(v1,W1,1%2,W2;A) is the pair distribution function in a 
system which the particles interact via a pairwise poten- 
tial 


un(r, w1,w2) = h”: (r) + dues (r, wi, we). (9) 


The density of Heisenberg particle py(r,w) can be split 
into the number density py(r) and a normalized factor 
On(r,w) as 


pu(t,w) = px(t)On(r,w) (10) 


where 


[t0ntr) =1 (11) 


In the local density approximation the Helmholtz free 
energy of HS mixtures system 


F™ = 5 | dedepu(e,.) [In(A},4rrpiz(r, w)) — 1112) 


j 5 I droy (r) [In(A%4mrpy(r)) — 1] 


rap oa 
Z ne) 

where 3 = 1/kpT is the inverse temperature, Aa is the 

thermal de Broglie wavelength of species a, the total den- 

sity p(r) = pp (r) + pn (r) and 7(r) = (7/6)p(r)o° is the 

packing fraction. 

In the mean field approximation where the pair 
distribution function takes its large-distance limit 
g(v1,01,%2,wWe;A) = 1 a2, the part of the free energy 
related to the spin-spin interactions in Eq. (0) becomes 


Fir = = zJ drıdrədwıdw2pH (r1, wy)ue> (r, W1, w2)PH (ra, wof)l3) 


After the decomposition py(r,w) = py(r)On(r,w), 
the minimum condition of the functional 
Qlox(r,w),pn(r)] is equivalent to the simultaneous 


minimization of the grand canonical functional with 
respect to the number densities 


dQ[px(r, w), pn(r), T, p] 


=0, 14 
Spa) = 
dQ(ox(r, w), pn(r),T, 4] =0 
dpn(r) 
and the orientational configuration 
b2 T 
[ex (r,w), pn(v), Hl = 0. (15) 
dOn(r, w) 
In the absence of the external field (VZ,(r,w) = 
VA,(r) = 0) the system is homogeneous in position, 


but it could be ordered in orientation. So we have 
pu(t) = pu, pn(t) = pw and Op(r,w) = On(w). In 
the mean field approximation we obtain the bulk expres- 
sion of the grand-canonical free energy density 


1 
y mr [On (v), px, PĒ: T, u] = f65(pH, PN) 


2 


here the Helmholtz free energy density of hard-sphere 
zi f23(pH, pn) is given by Carnahan and Starling 


Föslom pw) = [Inet 1] (17) 
An — 3n? 
peN [n(on AX) — 1] + - a = 7 
and 
Ss = [0n) s(w), (18) 
oS i dr4nr? J(r) = 8ne0°. (19) 


From the equilibrium condition and the mean field 
grand-canonical free energy we can obtain the equi- 
librium orientational distribution function 


eben Iss 


On) = Tipa Tes 


(20) 


Here we restrict to the case On (w) = Ow (0) and have 


1 
Ox(w) = On (0) = Lager acne. (21) 
AT 
where 
k(x) = In(z/ sinh z), (22) 


z = Boyd, (23) 


1. 
— sophia? + 22 | devon (w)inldrOn()] — ane — unon, 


with the average magnetization € = f dwOy (0) cos de- 
termined by 


1 


Then the grand function free energy density in an equi- 
librium state is expressed by 


E = coth(SpyJ'™€) — (24) 


1 
piimrlOn(8), px, en, T, u] = £65(eH, pw)(25) 
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IV. THE SPINODAL OF HEISENBERG HARD 
MIXTURE LIQUIDS 
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FIG. 1: (Color Online) A schematic plot of the eigenvec- 
tor which characterizes the phase transition. In the space 
(6p/p,6C,6€) the eigenvector can be described by angles 
@ and 0 which is in the range @ € (—90°,90°] and ¢ € 
(—180°, 180°]. When @ = 0 the phase transition is a pure 
ferromagnetic phase transition. The phase transition is pure 
condensation when @ = 90° and ¢ = 0. The phase transition 
is pure demixing when 6 = 90° and ¢ = 90°. 


At a stable equilibrium state the grand potential has 
its minimum and its variation with respect to the changes 
of number densities and magnetization should be positive 
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50/V = 


OpHOPN > 0. 


The variation of the grand potential can be rewritten 


in a matrix form 


1 
50/V = 5 ( PH Opn ôE ) MNH Mnn Myg Opn 
Men Men Mee ôt 
(27) 


We want to understand the character of the phase tran- 
sition, which combinations of condensation, phase sepa- 
ration, and ferromagnetic order fluctuations are leading 
to the phase transition. For the following analysis it is 
more convenient to rewrite the matrix in terms of the 
total density p = pp + py , the concentration C = py /p, 
and the ferromagnetic ordering, é. We reexpress the vari- 
ation of the grand potential in new fluctuations space 


(dp, 6C, 8E) 


Mpp Mpc Moe 


5n/V = 5 ( 6o/p oC ôE ) Mcp Moc Moeg ôC 
Mep Mec Mee ÔÈ 
(28) 


where total one particle density fluctuations dp = 
pu + 6pn, the demixing fluctuations 6C = ðL = 
Pp *|pndpy — puôpyn]. And the matrix M is defined by 


Mpp Mpc Mpg 
Mcp Mcc Mce (29) 
Mep Mec Meg 


M = 


which has three eigenvalues 1, A2,A3 with the corre- 
sponding eigenvectors X1, x2,x3. The positive eigenval- 
ues 1, Ag, and A3 of a stable state guarantee that the 
grand free energy will increase with respect to any varia- 
tion of total number density, magnetization, and concen- 
tration. If the smallest eigenvalue As vanishes, the system 
can deviate away from the original state without any in- 
crease of the grand free energy and becomes unstable, 
which is a phase transition. The eigenvector x, corre- 
sponding to the zero_eigenvalue characterizes this phase 
transition precisely ie and x, is the order parame- 
ter. Therefore, We investigate the direction of eigenvec- 
tor x, in the (69,6C,6€) fluctuations space, which can 
characterize the type of phase transitions. 

In Fig[] we show a schematic plot of the Euclidean 
eigenvector x, in the space (p/p, 8C, ÔE) with unit vec- 
tors (ez, €y,ez) of the Cartesian system. We normalize 
the eigenvector |x,| = 1 of the zero eigenvalue, and calcu- 
late the angles 0 and ¢ of x, in the spherical coordinate 


0 = arccos(x, : ez), (30) 


steel 3 
(=G =) (31) 


The angle ¢ and 0 describe the portion of the total 
number density, the concentration and the magnetiza- 
tion in the phase transition and is defined in the range 
—90° < 6 < 90°, —180° < ¢ < 180°. In general, the 
phase transition of the Heisenberg fluid is a combination 


ġ = arccos ( 


Muy Mun Mpe OPH 


p/p 


of condensation and ferromagnetic phase transition. If 
0 = 0, we have a pure ferromagnetic phase transition. 

or 0 near zero and ¢ ¥ 0, we have a ferromagnetic 

ominant phase transition accompanied by a weak con- 
densation and demixing fluctuations. If 6 = 90° and 
o = 0, we have then a pure gas-liquid phase transition. 
For 6 = 90° and ¢ near zero, we have a condensation 
dominant phase transition accompanied by a weak phase 
separation and ferromagnetic phase transition. If 0 = 90° 
and @ = 90°, we have then a pure demixing phase tran- 
sition. For 0 = 90° and ¢ near 90°, we have a demix- 
ing dominant phase transition accompanied by a weak 
ordering and condensation phase transition. Now we in- 
vestigate the phase behavior of the binary fluid by deter- 
mining the border of stable region from the zero point of 
the smallest eigenvalue x,. What we will obtain is the 
so-called spinodal. 


V. RESULTS AND DISCUSSION 
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FIG. 2: (Color Online) Plot of the eigenvalues 41, A2,and 
A3, vs T for the Heisenberg Hard sphere mixture fluid at the 
density p* = 0.8 and the concentration C = 0.8. The Curie 
point and the mixed spinodal point here are on the Curier 
line and the mixed spinodal of Fig. [B{f). 


In Fig[2] the eigenvalues A1, Az and A3 are shown as a 
function of the reduced temperature T* = 1/(e for the 
reduced density p* = po = 0.8 and the concentration 
C = 0.8. The eigenvalues Ag and A3 always keep to be 
positive, but the the smallest eigenvalue A, = A, ap- 
proaches zero when decreasing or increasing the temper- 
ature. The instability at higher temperature is actually 
on the Curie line where 6 = 0. And the smallest eigen- 
value at lower temperature meet the mixed spinodal line 
where 0 = 62° and ¢ = 80°. So we have a dominant 
phase separation accompanied by a weaker ordering and 
weakest gas-liquid phase transition. With the zero points 
of A, for different concentrations C we can get the spin- 
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FIG. 3: (Color Online) (a,d) The spinodal curve of the Heisen- 
berg Hard Sphere mixture fluid in different density sections 
(p* = 0.3,0.8). The blue dashed line, the Curies Line; the 
green solid line, the mixed spinodal curve; CEP, critical end 
point; IG, Isotropic Gas; FL, Ferromagnetic Liquid. (b,e) The 
magnetization € along the mixed spinodal curve with respect 
to concentrations. The pink solid line, €. (c,f) The angles 
along the mixed spinodal curve with respect to concentra- 
tions. The black solid line, the angle 0 ; the red dashed line, 
the angle ¢ 


odal phase diagram of the Heisenberg Hard sphere liquid, 
which is shown in Figi] 


A. Case: In the isotropic phases € = 0 


At the isotropic phase of the Heisenberg fluid, there is 
no total magnetization and x = 0. In this case we have 


_ O*k(x) 1 
oes (32) 
For all densities p the element Mj; is positive expect 
the element Mee, which can become zero and 6 = 0 cor- 
responds to a pure ferromagnetic phase transition. 
From Meeg = 0 we can get the Curie line of the Heisen- 
berg Hard mixture fluids. With the reduced density 
p* = po and the reduced temperature T* = 1/Ge we 
can express the Curie line as 
87 
=p" (33) 
3 
which agrees with the result of Li et al because of 
the simple mean field approximation. And the mean- 
field approximation that makes the phase transition be 


pure ferromagnetic. In a more accurate theory the phase 
transition is not pure ferromagnetic and should be ac- 
companied by the weak condensation and demixing phase 
transition. 


B. Case: € 40, the total density is fixed 


When € Æ 0 the determinant of the coefficient matrix 
M 


Mpp Mpc Moe 
Det Mcp Mcc Mce =0. (34) 
Mep Meo Mee 


which can be solved numerically. In Figi] (a), we show 
the Curie line (Blue Dash Line) and the mixed spinodal 
(Green Solid Line) of the HHNH mixtures with the fixed 
density p* = 0.3 section. 

the magnetization and the angle as a function of con- 
centrations. The phase diagram Will be discussed in de- 
tail and characterized the type of phase transitions in the 
method|19}. 

In Fig] (a) we obtain the spinodal phase diagram of 
the Heisenberg Hard spheres mixture liquid in the con- 
stant density in the temperature concentration plane. 
The mixed spinodal meets the Curie line at the criti- 
cal endpoint with the concentration C = 0.85 and the 
reduced temperature T* = 2.13. Below the temperature 
of the critical endpoint there is a first-order phase tran- 
sition between isotropic vapor and ferromagnetic liquid. 
The magnetization € along the mixed spinodal is shown in 
Fig[3](b) by using Eq.(24). The more Heisenberg particles 
the the liquid have, it show the bigger the magnetization. 
For the mixed spinodal the phase instability is a combina- 
tion of condensation, demixing and ferromagnetic phase 
transition. To characterize the phase instability precisely 
we investigate the angles ¢ and @ of the eigenvector x; 
along it shown in Fig] (c). For concentration near 1.0 
the angle ¢ is very small and positive and 0 > 45°, where 
the phase instability is predominantly related to the con- 
densation. The positivity of @ means that an increase of 
the total number density will be accompanied by a small 
demixing and ordering. This result is plausible, because 
the increase of concentration (more Heisenberg particles) 
will enhances the average attraction of the system which 
results in more ordering and demixing. The magnetiza- 
tion € for the concentration larger than 0.95 is larger than 
0.5. With the decrease of the concentration the angle ¢ 
increases but 0 decreases continuously. When the system 
approaches the critical endpoint C = 0.85, the angle 0 
approaches 0. The phase transition at critical endpoint 
is then a pure ferromagnetic phase transition in the mean 
field approximation. 

The spinodal curve and the CEP with C = 0.69 and 
T* = 2.89 at p* = 0.8 section is displayed in Fig[3] (d). 
Corresponding the magnetization in Fig[3](e) and the an- 
gles in Fig[3](f) have a obvious region 0.735 < C < 0.965 


that 0 > 45° and ¢ > 45°, in which the phase instability 
is predominantly related to the demixing phase transi- 
tion accompanied by the ordering and condensation. By 
comparing Fig[3] (c) with (f) the demixing predominant 
region disappear in the lower total density section. 


C. Case: #0 and Temperature is fixed 


0.70 90 


a b 
c 065 ~ FL 
k=] Ri 60 
3 CEP © 
5 0.60 Se g 
C IG Si D 
8 =s 2 30 J 
~ fa) 
5 0.55 ey — o 
E T=4.5 T=4.5 — o 
0.50 
0.80 0.85 0.90 0.95 1.00 0.80 085 0.90 0.95 1.00 
Density Density 
0.9 90 
ost \FL_-——_G d 
c o7 \, CEP 
2 0.6 IG ^, © 
m 0. 
= g 
zos ? 
2 04 A 3 
[e] om 
O 03 Ses 
on T=2.5 a 
03 04 05 06 07 08 09 1.0 03 04 0.5 06 07 08 09 1.0 


Density Density 


FIG. 4: (Color Online) (a,c) The spinodal curve of the Heisen- 
berg Hard Sphere mixture fluid in different Temperature sec- 
tions (T = 4.5, 2.5). The blue dashed line, the Curies Line; 
the green solid line, the mixed spinodal curve; CEP, criti- 
cal end point; IG, Isotropic Gas; FL, Ferromagnetic Liquid. 
(b,d) The angles along the mixed spinodal curve with respect 
to the total density. The black solid line, the angle @ ; the red 
dashed line, the angle ¢. 


The spinodal diagrams in the density concentration 
plane are shown in Fig. [4{a,c). In order to elucidate 
the differences of the mixture phase behavior at different 
temperatures, phase behavior is characterized by 0 and @. 
Fig. Ab) show in the T = 4.5 section the angles ¢ > 85° 
and @ < 45° along the mixed spinodal curve, in which the 
phase instability is predominantly related to the ordering 
phase transition accompanied by the weak demixing and 
tiny condensation. But in the lower temperature section 
T = 2.5 ( Fig[4{d) ) the instability predominantly related 
to the demixing accompanied by the weak ordering and 
gas-liquid phase transition at px > 0.6. Although the 
different temperature sections are quite similar topology, 
the predominant phase instability is considerable discrep- 
ancy between them. Recently, Fenzel at. presented con- 
centration pressure phase diagrams of the ideal Ising mix- 
ture fluids at different temperatures via the Gibbs ensem- 
ble Monte Carlo simulation [20}, which the phase behav- 
ior can qualitatively agree with the T sections calculated 
in our model. On the other hand the demixing phase 
transition in the mixtures can be easily observed in the 
lower temperature section but ferromagnetic phase tran- 


sition predominate in the higher temperature section. 


D. Case: € #0 and Concentration is fixed 
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FIG. 5: (Color Online) (a,d,g) The spinodal curve of the 
Heisenberg Hard Sphere mixture fluid in different Concentra- 
tion sections (C = 0.8, 0.9,0.99). The blue dashed line, the 
Curies Line; the green solid line, the mixed spinodal curve; 
CEP, critical end point; IG, Isotropic Gas; FL, Ferromagnetic 
Liquid. (b,e,h) The magnetization € along the mixed spinodal 
curve with respect to the total density. The pink solid line, €. 
(c,f,i) The angles along the mixed spinodal curve with respect 
to the total density. The black solid line, the angle 0 ; the red 
dashed line, the angle ¢. 


Density-temperature spinodal diagrams are plotted in 
Fig.()(a,d,g) where the different sections belong to dif- 
ferent values of the parameter C. In Fig[5] (a) the crit- 
ical endpoint is at p* = 0.345 and T* = 2.28 and the 
slope of the mixed spinodal is positive which was found in 
asymmetric binary dipolar mixtures|2]}. In Fig] (c) the 
phase instability is predominantly related to the demix- 
ing in the region p* > 0.55 and the ordering in the range 
of 0.345 < p* < 0.55. But the condensation fluctua- 
tions is always weakest in the C = 0.8 section. The 
positive slope of the mixed spinodal means that increase 
of the total number density, at fixed concentration, the 
first order phase transition temperature increase with in- 
creasing p*. When considering mixtures at C = 0.9 one 
finds from Fig. (5)(d) that the slope of the mixed spin- 
odal move towards flat. Since increase of the Heisenberg 
particle concentration implies that the condensation in- 
stability gradually predominate the phase transition, as 
shown Fig. (5)(f). In the limit C — 1.0 one recovers the 
MF spinodal diagram of the pure Heisenberg fluid con- 
tain only two fluid phases, which is an isotropic vapor 
and a ferromagnetic liquid [17]. Here we set C = 0.99, 
and then the spinodal diagram is shown in Fig[5](g) and 
the angle ¢ is always near the zero and @ is large than 45° 


when p* > 0.25 in Fig[5] (i). The predominant instabil- 
ity is related to the condensation with small ordering and 
demixing fluctuations which coincide with results of pure 
Heisenberg fluid. In particular, the angle is obtained as 
a function of the concentration in Fig] thereby provid- 
ing an estimate of the crossover between the demixing 
predominant and the condensation predominant phase 
transition as expected. Such the crossover phenomena 
can be precisely characterized by the method 


VI. CONCLUSIONS 


In this paper we have investigated the phase transition 
of HHNH mixtures by using the density functional theory 
in the mean-field approximation. The phase instability 
of the system is discussed with the method developed in 
Ref (19). From the matrix of the second derivatives of 
the grand canonical free energy Q with respect to the 
total particle density, the concentration and the magne- 
tization, we can determine the thermodynamic stability 
of the system. When the smallest eigenvalue of the ma- 
trix becomes zero, the system becomes unstable. The 
eigenvector corresponding the zero eigenvalue can char- 
acterize the phase instability precisely. In the total den- 
sity, concentration, magnetization space (dp/p, 6C, ôE) 
the normalized eigenvector can be described by two an- 
gles 0 and ¢ which is in the range —90° < 6 < 90° and 
—180° < ¢ < 180°, respectively. 

For temperature above of the critical endpoints, the 
angle 0 = 0 and the phase transition is pure ferromag- 
netic. This is the result in the mean-field approximation. 
With a more accurate approximation the angle 0 does 
not equal exactly to 0 and is just near 0. Then the total 
particle density and concentration will be related to the 
phase transition, which is predominantly ferromagnetic 
and accompanied by a weak condensation and demixing. 
Below CEP there is a first-order phase transition between 
isotropic gas and ferromagnetic liquid. The phase insta- 
bility along the spinodal near the ferromagnetic liquid is 
a combination of condensation, demixing and ferromag- 
netic phase transition. The phase behavior can topologi- 
cally be worthy to be compared with the ee ete on 
in the different T sections. Although the different tem- 
perature sections are similar topology, in the lower tem- 
perature section T = 2.5 the instability predominantly is 
related to the demixing but in the T = 4.5 section the or- 
dering phase transition is predominant. The slope of the 
mixed spinodal is the positive in the less concentration 
C = 0.8 sections and is negative in the C = 0.99 section. 
Hence by accurately investigating the angles we quanti- 
tatively describe that the phase instability continuously 
changes from predominant demixing phase transition to 
predominant gas-liquid phase transition. Although these 
crossover phenomena in the concentration sections may 
be considered as different phase transitions, our analysis 
shows that their origin in fact is same. 

These crossover phenomena within mean field approx- 


imation should be corroborated Monte Carlo simulations 
and further explore the combinations of phase transitions 
on the Curie line with more accurate approximation. 
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VII. APPENDIX: THE ELEMENTS OF THE 
MATRIX M 


In this appendix we give the expressions of the matrix 
elements (29) in detail, 


Mon = P x (C?Min +2CMan (35) 


— 2° Myn + (1 — C)’Mww), 


Moc = #(Myu —2Myvn + Myyn), 36) 
Mpg = Mep = pCMue, 37) 

Mog = Mec = pMue, 38) 

Mop = Mpc = 0? x (CMun 39) 


+ (1-2C)Muy ~(1-C)Myy), 


where the elements of the matrix are related to 


Ək 
Muu = BJ? intP pul” ðr a 2 + 2Jinté one aa) 
2 p 4n-3n? 
pl ee E 
BpH opt i 
1 PE n 
M = — p bUm’ 41 
NN HoN + Doron (41) 
k(x 
Mee = BIinePiz xt + Jint Pir, (42) 
AS 
1—n 
M = Myny = 43 
HN NH Dawe (43) 
3? k(x) 
Mue = Men = bJ PHE Anda (44) 
Ok 
+ Ma 2 int PHE, 
Myne = Men = 0. (45) 
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